Signature of Sub GeV Dark Matter particles at LHC and TEVATRON 
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In this letter, we investigate the production of light dark matter particles at LHC in light of the model (N = 
2 SUSY inspired) proposed in Ref. (l[ and demonstrate that they will be copiously produced if the colored 
messengers F q are lighter than 1 TeV. We expect up to 10 6 events if ~ 500 GeV, assuming a ~1 fb -1 
luminosity. In addition, we show that, even if > O (1) TeV, searches for F q production at LHC are promising 
because a kinematical signature can be used to separate the signal from background. This signature is similar 
to that expected in supersymmetric scenarios. Hence, our study shows that most of the mp range could be 
constrained using LHC data. This should encourage further studies since they could infirm/confirm the MeV 
DM scenario. 



I. INTRODUCTION 

It is often argued that an important motivation for carrying 
new physics searches at LHC and TEVATRON is the possibil- 
ity of discovering new particles that could enlight the nature of 
dark matter (DM). If DM is made of thermal particles, its rel- 
atively small abundance today then guarantees a coupling to 
Standard Model particles and therefore suggests that it could 
be produced in high energy experiments. 
Unfortunately DM indirect detection at LHC and TEVATRON 
is not easy. The DM signatures that have been studied for 
a broad range of models (see for example Supersymmetry, 
lepto-quarks, Kaluza-Klein, little Higgs) involve e.g. the pro- 
duction of leptons, di-leptons, jets as well as missing energy 
J2]. In these frameworks, the DM mass is generally assumed 
to lie from a few GeV to a few TeV. Yet this is only a sub- 
set of all the DM possibilities that have been proposed in the 
literature. Besides, there is no evidence for such GeV- TeV 
DM particles in laboratory experiments so far (apart, perhaps, 
from JH01), despite intensive searches. In fact, direct detec- 
tion experiments have now demonstrated that, if heavy neutral 
particles exist with a mass in the 100 GeV- TeV range, their 
interactions with ordinary matter must be small, e.g. 
implying fine-tuned solutions (e.g. 10]) to simultaneously ex- 
plain the DM relic density. 

In this context, it is worth investigating other DM candi- 
dates. Here we shall concentrate on 1-100 MeV particles 
lUl H which were first introduced while studying the effect of 
DM interactions (whatever the nature of DM) on Large Scale 
Structure formation JUGi!. A subclass of these candidates 
(DM with a mass of a few MeV) has received a lot of atten- 
tion after it was realized that it could explain the emission of 
511 ke V photons in the galactic centre while conventional as- 
trophysical sources generally fail to reproduce the observation 
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The phenomenology and detectability of such MeV particles 
have been studied in detail 01 [l2Hl4ll . However, all previ- 
ous studies assumed low energy experiments (where chances 



of detection seemed greater owing to the very small DM cou- 
pling to electrons 11511 ). Here we show that similar studies for 
high energy experiments are actually much more promising 
(on the reasonable assumption that DM is coupled to quarks). 
The paper is organized as follows. In Sec. [II] we introduce 
our model. In Sec. [HI] we compute the number of events ex- 
pected in an experiment with one fb^ 1 luminosity. Motivated 
by previous analysis of similar final states, we study the kine- 
matical characteristics associated to our signal and compare 
it with those of the background in Sec. QV] We conclude in 
Sec.E 



II. THE MODEL 

Motivated by the model proposed in Ref. {!]] (based on N = 
2 supersymmetry), we consider a Lagrangian L containing a 
term 
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where 8, 7 - denotes the Flavour indices, S an SU(2) singlet 
scalar field (our DM candidate), F a heavy fermion and / a 
Standard Model (SM) particle. 

Such a term looks like the "standard" coupling between a 
sfermion, a SM fermion and a neutralino in N = 1 supersym- 
metry, except that in our model F is a spin 1/2 particle (asso- 
ciated with f) and S is a neutral scalar. To preserve gauge in- 
variance, Fr is a doublet of SU (2) while Ft is a singlet fll [l6ll . 
Hence terms like SFlFr or HFtfg (with H the Higgs field) are 
absent from the Lagrangian since they are not gauge invari- 
ant. Owing to F quantum numbers, one could introduce other 
terms, such as HFlFr or HFtfi- The former is expected to 
give a mass to the F (hereafter denoted by mp). However, 
the exact value of mp may arise from the existence of several 
Higgs doublet in the underlying theory and, potentially, con- 
tributions from soft symmetry breaking (like inN ~l SUSY). 
The later term introduces a mixing mass matrix between the 
SM fermion and the F. However we assume the existence of 
a new symmetry (M-parity) Hl flrjll which kills such a term. In 
principle, a phenomenological study ofN = 2 supersymmetry 
would require to write the full Lagrangian and would lead to 
all difficulties encountered already in N = I SUSY. However 
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we argue that, with this term only, one can already learn about 
such theories. As an illustration, we focus on the MeV DM 
scenario (which was based on this Lagrangian JH [HI [l5lo . 
However a similar analysis could be done for heavier DM par- 
ticles. 

Let us now focus on the quark sector. For simplicity, we as- 
sume that all F q have a mass mp = mp ^ m q and all the F q 
couplings to their corresponding SM quark have the same val- 
ues. This may cancel potential contributions to rare meson de- 
cays. In addition, we do not introduce any flavour mixing be- 
tween the various F. CP phases are set to zero. Owing to these 
properties, there should not be any large FCNC contribution 
in this set up. We make similar assumptions for the leptonic 
sector, except that ^ mp v (y^ mp), DM pair annihilation 
into neutrinos could thus insure the correct relic density lfl7tl 
(even though the relic density criteria may be alleviated by the 
assumption of two DM particles iflHl ). 

We can now study the signatures associated with this scenario. 
To illustrate our purpose, we consider couplings c/ , varying 
between [0.3,3]. Couplings above unity may, of course, ap- 
pear rather unlikely as they may induce unseen anomalies in 
particle physics measurements (depending of mp • ) and pro- 
duce a very bright, yet extremely narrow, monochromatic line 
through quark box diagrams at an energy E = m.d m in our 
galaxy lfl9tl . However, for this analysis, we use these very 
large couplings as benchmark points to determine the typical 
values of the MeV DM coupling to quarks that can be probed 
at LHC. Note that the F q being "colored", they can be directly 
produced through gluon-gluon fusion. 

III. F q PRODUCTION CROSS SECTIONS IN 
PROTON-PROTON COLLISIONS 

A. F q production in qq, qq, gq collisions and gg fusion 

The qq — > F q F q cross section involves a t and u-channel DM 
exchange as well as a s-channel gluon exchange. The qq — > 
F q F q channel is similar but there is no s-channel gluon ex- 
change i20ll . The qg — > F q S process involves a t-channel F q ex- 
change and a quark s-channel exchange while the gg — > F q F q 
is based on gluon fusion and t + u-channel F q exchange. The 
behaviour of all these cross sections with respect to the parton 
energy is displayed in Fig. Q] (left panel). 
Interestingly enough, the F q production cross section through 
qq collisions remain significant even for mp q > 300 GeV and 
ci ,- < 3. For example, it is about a few 30 pb for mp = 300 
GeV, ci r = 1 and y/s~2 mp q . This can be understood by com- 
puting the matrix squared amplitude and writing the centre of 
mass energy as s = 4 m 2 F (1 + e 2 ) (where e can be very large). 
For very large or very small values of e, one finds that o=< 1 /e. 
But the cross section has a maximum when e satisfies: 

-4 [7(4 + c 4 r ) + 12c 2 c 2 ]£ 6 - 16 [3 (c? + 4) + 5 c 2 c 2 ]e 4 

-3 [5 (cj + c?)+8c 2 c 2 ]e 2 +2 [3(4 + 4)+5c 2 c 2 ]=0 

(assuming cos 8 = 0). Thus, if c/ = c r = 1, O qq -^p p is max- 
imal at e « 0.47 (or e « 0.64, if one integrates over cos 8 in- 
stead of setting it to zero). We then find a « 30pb, if mp = 300 



GeV. This cross section is therefore very large when the F q can 
be produced almost on-shell. 



B. F q production cross sections in p-p collisions 

To compute these cross sections, we used the function hCol- 
lider implemented in micrOMEGAs l2lll . We have checked 
that our results were consistent with our analytical expressions 
(before convoluting with the parton distribution function) and 
the output from the VEGAS and Easy 2x2 from CalcHEP 
JH. Results from VEGAS and hCollider differ by a factor 
~ 1 .5 at high energy but this is due to the fact that we did not 
into take into account some QCD corrections in VEGAS. 
Results are displayed in Fig. Q] (middle and right panel). For 
not too large c/ ,- couplings, the dominant cross section corre- 
sponds to the pp — > F q Fq process since it involves gluon fu- 
sion and gluon exchange process. The latter can reach up to 
a few nb when s/s > 9 TeV (assuming nip = 300 GeV) but it 
rapidly falls off with mp q . For y/s = 10 TeV and mp = 2 TeV, 
it is only about 10 pb (see right panel of Fig . |T} . 
These numbers suggest that TEVATRON may set a limit on 
mp . Due to interference between the S and gluon exchange 
diagrams, the pp — > F q F q production cross section is maxi- 
mal for C{ = c r = 3 and minimal for C[ = c r ~ [1,2]. For 
ci = c r < 1, the gluon exchange is the dominant process. 
Hence, the cross section is fixed by QCD couplings. As a 
result, we can set a limit on the only free parameter that is 
left, i.e. mp q . By analogy with LeptoQuark (jets+neutrinos) 
searches, we found that mp should be greater than 450 GeV, 
see Figj3] However, this limit may also be 300 GeV, given 
that searches for this very model have not been implemented 
in TEVATRON experiments yet. Hence, to be "conservative", 
we shall use mp q = 300 GeV in the next section. Any greater 
value of mp q will imply a smaller cross section. 

IV. EVENT SIGNATURE AND BACKGROUND 

As shown in Sec. Hill F q should be produced significantly in 
proton colliders at high energy. In Fig. |2] (left) we present 
the production cross sections as a function of ci,c r parame- 
ters. The production setup is pp collisions at y/s = 7 TeV. 
These are expected to be the initial LHC conditions for the 
beginning of the physics program. Our goal is to study the 
kinematical properties of the signal and compare them to the 
associated background in order to establish the discovery po- 
tential of this signature at the LHC, particularly with the AT- 
LAS detector 112311 . ATLAS will be able to measure different 
observables on the objects that compound the final state of our 
signal, namely: jet identification and missing energy. Here a 
jet is understood as the imprint left by the hadronization pro- 
cess of a high energy quark in the detector material. 
Since each F q decays into a jet + missing energy {fir), the 
associated background would be composed by: First, Z+jets 
where Z decays into two neutrinos. Second, tt where t decays 
into W, b and W into /,v. In cases where the lepton falls out 
of the region of possible identification (r| > 2.8 for ATLAS), 
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F„ production cross sections at parton level 



w 0.4 
c 

f 0.35 
w 

$ 0.3 

o 

°0.25 
0.2 
0.15 
0.1 
0.05 




uu-» FF 
uu — » FF 



98-»FF 



j D. Albomoz^asquez, C Boehm. J. Idarraga j 



2 3 4 5 

Parton momentum [TeV] 



F u production cross sections in p-p collisions 
£■ 3.5p 



F u prodjction cross sections in p-p collisions 



3 4 5 6 7 
Proton momentum [TeV] 




■ I I I I 1 1 1 1 1 1 


1 






PP -> FF 






— pp-»FS 






— pp-» FF 




r "v\ 
r "sjv 






! D. AlbOrnbz Vasq 


ez, c. Boehm, J. Iflaitfiga 1 


- 







1500 2000 2500 
m F in [GeV] 



FIG. 1: F q production cross sections. Left panel: at the parton level, assuming mp = 300 GeV, c; r = 1 and obtained by using the function cs22 
implemented in micrOMEGAs. Middle panel: at the proton level, with trip = 300 GeV, c\ r = 1 and obtained by using hCollider. Right panel: 
dependence of the cross section with mp . Same parameters as middle panel but ^/s = 10 TeV. 



this signature can mimic the signal, tt will be produced with 
an enormous cross section at the LHC (see Table for all cross 
sections) and it has to be taken into account in realistic sim- 
ulations. Even though, previous studies have shown that it is 
not a real competitor for signals with a final state composed 
by 2 jets + f.j as shown in |2S](p.l595), or other combinations 
of leptonic, hadronic plus fir final states as studied in l24ll . 
One of the advantages is that ATLAS has fb-tagging capabil- 
ity, and whenever an event clearly contains a /b-quarkjet, we 
can reject it. Third, the SM known processes involving WZ, 
ZZ and WW production. In the first and second case it can 
perfectly mimic the signal when W or Z decays into two jets 
and Z decays into neutrinos. It has been proved f2(|(p.l595) 
that a cut based analysis based on the kinematical properties 
of signal an backgrounds can eliminate this background com- 
pared to signals down to a few fb. The same techniques are 
successful against the Z+ jets background mentioned above. 
WW, like in the case of tt, can only reproduce the signal when 
one lepton falls out of the detectable region. The cross sec- 
tions in the table below are production cross sections, and do 
not include the branching ratio. These cross sections estima- 
tions have been produced with with the MadGraph genera- 
tor l25ll including pre-selection cuts compatible with ATLAS 
calorimetry and tracking, according to those used by the AT- 
LAS collaboration in 12J] . 

Even in pessimistic scenarios, where smaller values of the 
dark matter couplings to quark (for example C[ = c r ~ 0.3) 
lead to cross sections of the order of a few pb for mp < 200 
GeV (or down to 40 fb for mp < 500 GeV), the signal identi- 
fication with this type of final state is possible as shown in J2l- 
In Figj2] (right), we take the most competitive background, 
Z+jets (in this case we consider only Z + 2 jets at parton level 
produced with the MadGraph generator [25]) and plot the pj 
and r| distributions against the same observables associated to 
one of the F q s in the signal. The hardness of the pj distri- 
bution for the signal, as opposed to that of the shown back- 
ground, is a typical characteristic of a two body decaying sig- 
natures at high energy, and we include it here to show that we 
have successfully implemented a MonteCarlo (MC) machine 
which allows us to study further the kinematics of this signal. 



We are going into a full simulation of the associated final state 
using the Geant-4 simulation of the ATLAS detector already 
thoroughly tested by the ATLAS collaboration JH] 0]. 
Comparing the kinematical characteristics of our signal 
(which we have been able to study for the first time) with 
the signature of different signals which share the same final 
state (jets+missing energy) and are known as candidates for 
discovery at the LHC, we claim that current and forthcoming 
high energy experiments should have the ability to constrain 
the scenario presented in this letter. 
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V. CONCLUSION 

In this letter, we investigate the case of DM particles coupled 
to SM quarks through a new type of colored charged fermions 
Fg ( as predicted in models yj] inspired from N = 2 SUSY 
iflrjll ) and study F q production at LHC and TEVATRON. We 
show that, in the case of sub GeV DM, up to 10 6 events 
could be produced in a collider with a 1 fb~ 1 luminosity if 
mp q ~ 300 GeV, and about 10 3 events if mp ~ 2 TeV (with 
t/s =10 TeV). In addition we found that the kinematic char- 
acteristics of the signal in contrast to those of the associated 
backgrounds can be used on a based cut analysis on simulated 
or real data for mFq > O(l). This indicates that searches for 
this type of couplings (at least in the case of light DM parti- 
cles) should be possible at LHC and TEVATRON and should 
definitely motivate further studies in high energy experiments. 
One should also remember that a large production of unstable 
colored F q particles could lead to a large production of muons 
(after hadronization of the jets in the detector), which could be 
very useful to constrain such a scenario. This will be inves- 
tigated in a forthcoming study. Finally, since the kinematical 
signature of this model is similar to that expected in some su- 
persymmetric scenario, a study of the spin of the F q particle 
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FIG. 2: c;,c> parameter's phase-space giving the production cross sections at the LHC at y^s = 7 TeV (left). Kinematic characteristics of the 
signal as opposed to those of the Z boson in the associated Z+jets background (right). In the figures (right), we have normalized the signal and 
background cross sections so that they become comparable. We use c/ = c r = 1, mp = 300 GeV, */s = 7 TeV. 



F u production cross sections at TEVATRON 
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FIG. 3: F q production cross sections at Tevatron. 



may be required to help for the identification of the dark mat- 
ter. 
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